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We propose a method to control the number of species of lattice fermions which 
yields new classes of minimally doubled lattice fermions. We show it is possible 
to control the number of species by handling 0{a) Wilson-term-like corrections in 
fermion actions, which we will term "Twisted-ordering Method" . Using this method 

we obtain new minimally doubled actions with one exact chiral symmetry and exact 
(N 

locality. We classify the known minimally doubled fermions into two types based on 
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the locations of the propagator poles in the Brillouin zone. 
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I. INTRODUCTION 

Since the dawn of lattice field theory [l|, the doubling problem of fermions has been one 
of the notorious obstacles to QCD simulations. The famous no-go theorem of Nielsen and 

n 

Ninomiya [2[ most clearly reveals the fate of lattice fermions. It states that lattice fermion 
actions with chiral symmetry, locality and other common features must produce degrees of 
freedom in multiples of two in a continuum limit. In contrast, phenomenologically there exist 
only three quarks with masses below the underlying scale of QCD. By now several fermion 
constructions to bypass the no-go theorem have been developed, although all of them have 
their individual shortcomings. For example, the explicit chiral symmetry breaking with 
the Wilson fermion approach [11] results in an additive mass renormalization, which in turn 
requires a fine-tuning of the mass parameter for QCD simulations. Domain-wall fermion 
[3|, |4j and overlap fermion p, |6[ do not possess exact locality, leading to rather expensive 
simulation algorithms. These approaches attempt to realize single fermionic degrees of 
freedom by breaking the requisite conditions for the no-go theorem. 

On the other hand, there is another direction to approach numerical simulations. Ac- 

n 

cording to [7j], Hypercubic symmetry and reflection positivity of actions result in 2 d species 
of fermions where d stands for the dimension. Thus it is potentially possible to reduce 
the number of species by breaking hypercubic symmetry properly. Actually, the staggered 



fermion approach [8Hl0l|. with only 4 species of fermions does this and possesses flavored- 
hypercubic symmetry instead. However this requires rooting procedures for the physical 2 
or (2 + l)-flavor QCD simulation, which is known to mutilate certain processes. 

An interesting goal is a fermion action with only 2 species, the minimal number required 
by the no-go theorem. More than 20 years ago, Karsten and Wilczek proposed such a 



minimally doubled action ll|, Il2j . Recently one of the authors of the present paper 
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15| has also proposed a two-parameter class of fermion actions, inspired by the relativistic 



condensed matter system, graphene 16j. These minimally doubled fermions all possess one 
exact chiral symmetry and exact locality. As such they should be faster for simulation, at 
least for two-flavor QCD, than other chirally symmetric lattice fermions. 
However it has been shown 
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20] that we need to fine-tune several parameters for a 



continuum limit with these actions. This is because they lack sufficient discrete symme- 



try to prohibit redundant operators from being generated through loop corrections 2ll-l24|. 



Thus the minimally doubled fermions have not been extensively used so far. Neverthe- 
less, there is the possibility to apply them to simulations if one can efficiently perform the 
necessary fine-tuning of parameters. To understand how far we can reduce the number of 
fine-tuning parameters, it should be useful to investigate and classify as many minimally 
doubled fermions as possible. 

In this paper we propose a systematic method to reduce the number of doublers, which 
we term the "twisted-ordering method". In this way we obtain new classes of minimally 
doubled fermions. We study symmetries of these fermion actions and show CP invariance 
and Z4 symmetry. It is also pointed out that they all also require fine-tuning of several 
parameters for a correct continuum limit because of lack of sufficient discrete symmetries. 
We classify the known minimally doubled actions into two types. By this classification we 
can derive several further minimally doubled actions deductively. 

In Sec. [Til we present the twisted-ordering method and derive a new class of minimally 
doubled fermion actions. In Sec. IIIII we show that two classes of minimally doubled actions 
arise from the original twisted-ordering actions. In Sec. IIVI we investigate the symmetries of 
these fermion actions. In Sec. |V]we discuss two classifications of minimally doubled fermions. 
Section I VI I is devoted to a summary and discussion. 

II. TWISTED-ORDERING METHOD 

In this section we propose a systematic way of controlling the number of species of lattice 
fermions within the requirement of Nielsen- Ninomiya's no-go theorem. We will firstly discuss 
the 2-dimensional case to give an intuitive understanding of this mechanism. Then we will 
go on to the 4-dimensional case and show we can construct a new minimally doubled action 
by this method. 

A. 2-dimensional case 

Let us begin with the following simple 2-dimensional Dirac operator in momentum space 
with 0(a) Wilson-like terms. 

D(p) = (sinpi + cospi — 1)271 

+ (sinp 2 + cosp 2 - 1)«72, (1) 
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FIG. 1: Zeros for the 2-dimensional untwisted case are depicted within the Brillouin zone. Red- 
dotted and blue-solid lines stand for zeros of the coefficients of 71 and 72 respectively in Eq. (fip. 
Black points stand for the four zeros of the Dirac operator. 

where px, P2 and 71, 72 stand for 2-dimensional momenta and Gamma matrices respec- 
tively. Here the deviation from the usual Wilson action is that the 0(a) terms are ac- 
companied by Gamma matrices. There are four zeros of this Dirac operator at (pi,p 2 ) = 
(0, 0), (0, 7r/2), (tt/2, 0), (tt/2, tt/2). This means that the number of the species is the same 
as for naive fermions. 

Next we "twist" the order of the 0(a) terms, or equivalently, permute cospi and cosp2 
to give 



D(p) = (smpi + cosp 2 — 1)271 
+ (sinp 2 + cospi - 1)272- 



(2) 



Here only two of the zeros ((0, 0) and (tt/2, it/ 2)) remain and the other two are eliminated. 
Thus the number of species becomes two, the minimal number required by the no-go theorem. 
As seen from this example, twisting of the order of 0(a) terms reduces the number of species. 
We call this method "twisted-ordering" in the rest of this paper. We depict the appearance 
of zeros in the Brillouin zone for the above two cases in Fig. [1] and [2j The red-dotted and 
blue-solid curves represent the zeros of the coefficients of 71 and 72 respectively. These 
curves intersect at only two points for the twisted case in Fig. [2j 

Now consider excitations and the physical (75) 's at the two zeros. Here the physical 75 
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FIG. 2: Zeros for the 2-dimensional twisted case is depicted within the Brillouin zone. Red-dotted 
and blue-solid curves stand for zeros of the coefficients of 71 and 72 respectively in Eq. ([2]). Black 
points stand for two zeros (0, 0) (vr/2, vr/2) of the Dirac operator . 

at a zero is defined as the product of four coefficient matrices of momenta in the excitation 
as we will concretely show later. Now the Dirac operator is expanded about the zeros as 



D {1 \q) = hiqi+ i72?2 + 0(q 2 ), 
D {2 \q) = -ijiq 2 - «72<?i + 0(q 2 ) 



(3) 
(4) 



Here we expand with respect to g M defined as p^ = p^ + g M and we denote the two expansions 
as D*- 1 ' for the zero (0,0) and D^ for the other at (•7r/2,7r/2). Here momentum bases at 
the zero (0, 0) are given by b-j/ = (1, 0) and b 2 = (0, 1) while those at the other (7r/2,7r/2) 
are given by b^ = (0, —1) and b 2 = (—1,0). This means excitations from the two zeros 
describe physical fermions on the orthogonal lattice. And the 75 at p = (0,0), which we 
denote as 75 , is given by 

(5) 



75 



(i) 



7i72, 



while the 75 at p = (vr/2, vr/2), which we denote as 7 5 , is given by 



(2) ; c 



(2) 
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(i) 



(6) 



This sign change between the species is a typical relation between two species with minimal 
doubling since Nielsen-Ninomiya's no-go theorem requires fermion pairs to possess chiral 
charges with opposite signs. 



B. 4-dimensional case 

Now let us go to the 4-dimensional twisted-ordering method. Here we again start with 
the following simple action with 0(a) Wilson-like terms 

D(p) = (smpi + cospi — 1)271 

+ (sinp 2 + cosp 2 - 1)«72 
+ (sinp 3 + cosp 3 - 1)273 
+(sinp 4 + cosp 4 - 1)274. (7) 

Here each of the four momentum components can take either or tt/2 for zeros of this Dirac 
operator. It means the associated fermion action gives 16 species doublers, the same as for 
naive fermions. 

Then we firstly consider the case with twisting the order of 0(a) terms once 

D(p) = (sinpi + cospi — 1) 271 

+ (smp 2 + cosp 2 - 1)272 

+ (sinp 3 + cos£> 4 - 1)273 

+ (sinp 4 + cosp 3 - 1)274, (8) 

where we permute the order of cosp3 and cosp 4 . Here we find that there are 8 zeros of the 
Dirac operator at 

p 1 = 0oi7i/2, p 2 = 0om/2, (p 3 ,p 4 ) = (0, 0) or (tt/2, tt/2). (9) 

The associated action yields 8 species of doublers. As seen from this example, the twisted- 
ordering method reduces the number of species also in 4-dimensions. In this case twisting 
is performed only once, and the number of species is reduced by a factor of two. 
Secondly we permute cospi and cos£>2 in addition to cos£>3 and cosp4 to give 

D(p) = (sinpi + cosp 2 - 1)271 

+ (sinp 2 + cospi - 1)272 

+ (sinp 3 + cosp 4 - 1)273 

+ (sinp 4 + cosp 3 - 1)274. (10) 



In this case there are only 4 zeros at 

(p u p 2 ) = (0,0) or (tt/2, tt/2) (p 3 ,p 4 ) = (0,0) or (tt/2, tt/2). (11) 

Thus the associated action yields 4 species of doublers. This example indicates that more 
twisting results in additional reduction of species. In these two examples, Eqs. (jHJ) and ffTUj) . 
twisting is performed partially, and so we call this kind of procedures "partially-twisted- 
ordering" . 

Finally we consider a case that the order of 0(a) terms are maximally twisted: We 
permute the cos terms in a cyclic way to give 

D(p) = (sinpi + cosp 2 — l)*7i 
+ (sinp 2 + cosp 3 - l)ij2 
+ (sinp 3 + cosp 4 - 1)273 
+ (sinp 4 + cospi — 1)274. (12) 

In this case there are only 2 zeros at 

(pi,P2,P3,P4) = (0,0,0,0), (vr/2,7r/2,7r/2,vr/2). (13) 

We have confirmed there are no other real zeros of this operator numerically, and we will 
discuss an analytical proof in the next section. There are only two species in this case, the 
minimal number required by the no-go theorem. This is a new type of minimally doubled 
fermion on the orthogonal lattice. 

Now we study the excitations and 75 's at the two zeros in a similar way to the 2- 
dimensional case, 

D {1 \q) = 27191 + 272^2 + 273<?3 + 274<?4 + 0(q 2 ), (14) 

D (2 \q) = -27i02 - «72<?3 -2'73<?4- 2'74<?i + 0(q 2 ), (15) 

where we expand the Dirac operator ( TT2|) around zeros with respect to g M defined as p^ = 
p^ + q/j, and we denote the two expansions as D^- 1 ' for the zero (0, 0, 0, 0) and D^ 2 ' for the other 
(7r/2, 7r/2, 7r/2, tt/2). Momentum bases at the zero (0, 0, 0, 0) are given by b^ = (1, 0, 0, 0), 
b^ 1} = (0,1,0,0), b 3 X) = (0,0,1,0) and b 4 1} = (0,0,0,1) while those at the other zero at 
(tt/2, vr/2, tt/2, tt/2) are given by b^ = (0,0,0,-1), bf = (-1,0,0,0), b 3 2) = (0,-1,0,0) 
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(2) 

and h>4 = (0, 0, —1,0). This means excitations from the two zeros describe physical fermions 
on the orthogonal lattice again. And the 75 at p = (0,0,0,0), which we denote as 75 , is 
given by 

75 = 71727374, (16) 

while the 75 at p = (vr/2, n/2, n/2, vr/2), which we denote as 75 is given by 

( 2 ) (!) /i7\ 

75 = 74717273 = -75 ■ ( 17 ) 

This is a typical relation between the two species with minimal-doubling. 

Gauging these theories is straightforward. One merely inserts the gauge fields as link 
operators in the hopping terms for the action in position space. Specifically, for the twisted- 
ordering minimally doubled fermion in position space we have 

n,fj, 

(18) 

where we define 

fl,2,3 (/x = 2,3,4) 
H-l=\ (19) 

{ 4 (/i=l). 

We note that in the twisted-ordering method, we control the number of species by choosing 
the extent of breaking of discrete symmetries of the action. Twisting of 0(a) terms leads to 
breaking of the hypercubic symmetry. On the other hand it is well-known that hypercubic 
symmetry of the lattice fermion action requires 16 species as shown in |7J. Thus it is 
reasonable that more breaking of discrete symmetries can lead to less doubling. 

III. MINIMALLY-DOUBLED FERMIONS 

In this section we discuss two classes of minimally doubled actions obtained from the 
original twisted-ordering action ( TT2l) . The first one, which we call the "dropped twisted- 
ordering action", is constructed as following: We drop one of cosp^ — 1 terms in the Dirac 



operator (TT21) . for example, drop cospi — 1 as 

D{p) = (siripx + cosp 2 — l)«7i 
+ (smp 2 + cosp 3 - 1)^72 
+ (sinp 3 + cosp 4 - 1)^73 
+ sinp 4 27 4 . (20) 

This has two zeros given by 

(pi,P2,Pa,Pi) = (0,0,0,0), (vr, 0,0,0). (21) 

Thus the associated action is also a minimally doubled action. Here the two (75) 's at the 
two zeros also have opposite signs 

7 5 X = 71727374, 

7? } = (-7i)727374 = -75 , (22) 

where 75 corresponds to (0, 0, 0, 0) and 75 to (7r, 0, 0, 0). We can also confirm this action 
describes physical fermions on the orthogonal lattice by looking into excitations from the 
zeros as in the previous section. We depict the appearance of zeros for the two dimensional 
case of the dropped twisted-ordering action in Fig. |3j We see there are only two zeros as 
(0,0) and (vr, 0). The 4-dimensional case is similar. 

Here we comment on definition of the equivalence of actions. We can convert one fermion 
action to an equivalent one with a different form by momentum shift and rotation of Gamma 
matrices. Here the momentum shift means 

p M ->• ±p ll + k h , (23) 

where k^ are some constant four vector. The rotation of Gamma matrices means 

7m ->■ C ^lu, (24) 

where C^ v is an orthogonal matrix. Although by using such a shift and rotation we can 
control the location of two zeros with the distance and direction between them fixed, the 
original and transformed actions are clearly equivalent. Now we can easily see the original 
( TT2l) and dropped twisted-ordering (1201 actions never translate into Borici [l4j or Karsten- 
Wilczek [111. Il2| actions through these two procedures. It means our actions are independent 
and new ones. 
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Dropped case 
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FIG. 3: The two zeros for the Dirac operator (|20p in the 2-dimensional analogy. Red-dotted and 
blue-solid lines stand for zeros of the coefficients of 71 and 72, namely (sinpi + cosp2 — 1) and sinp2 
respectively. Black dots show two zeros (0,0) and (vr,0). 

Next we consider a second generalization where we turn on a parameter in the action 
( fT2|) as following, 



D{p) 


= (sin pi + cosp 2 _ 


- a) z'7i 




+ (sin£>2 + cosp 3 - 


- a) ?72 




+ (sinp 3 + cosp4 - 


- «)«73 




+ (sinp 4 + cospi - 


- a) ^74 



(25) 

Here we replace unity with a positive parameter a in the constant terms of the operator. We 
will also call this one the twisted-ordering fermion action in the rest of this paper. The a = 1 
case corresponds to the original case. Now, in order to investigate the parameter range within 
which "minimal-doubling" is realized, we define sinpi = x, sinp2 = y, sinp^ = z, sinp^ = w 
and rewrite the equations for zeros in terms of these variables 



x 



± y/l - y 2 = a, 

y± Vl - z 2 = a, 
z ± Vl — w 2 = a, 



w 



±VT 



x* = a. 



(26) 
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To make it easier to see we rewrite this as 

x 2 + y 2 — 2ax = 1 — a 2 , 
y 2 + z 2 — 2ay = 1 — a 2 , 
z 2 + w 2 — 2az = 1 — a 2 , 
w 2 + x 2 - law = 1 - a 2 , (27) 

where in principle there are 16 complex solutions. The number of solutions of this system 
is the same as the original ones ( |25|) ff26|) . Thus we can study the range of the parameter for 
minimal- doubling by looking into how many real solutions this system of equations has. By 
eliminating y, z and w, we obtain the following equation for x 

a m (2x 2 - 2ax + a 2 - l) 8 = 0. (28) 

where we can confirm this by using software such as Maxima. This means, except for a = 0, 
there are only the two solutions 



a ± \/2 — a 2 
x = y = z = w= . (29) 

From this a should be below a = a/2 so that the solutions are real. For a = the 
above equation ( 128|) is satisfied by any value of x, which means that zeros appear as a one- 
dimensional curve, not a point. We will discuss this in detail later. We find there are two 
real solutions in the range < a < v2, which corresponds to a pair of minimally doubled 
species. This is the range of the parameter a for minimal-doubling. We can also verify this 
minimal-doubling range numerically. 

What is significant for a = is that zeros of the Dirac operator (1231) appear as an one- 
dimensional curve, not a point, in the momentum space. We can confirm this as following. 
Equations for zeros of the Dirac operator (125]) at a = are given by 

sinpx + cosp 2 = 0, 
sinp 2 + cosp 3 = 0, 
sinp 3 + cos£>4 = 0, 
sinp 4 + cospi = 0. (30) 
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These equations are satisfied by 

p 2 = ±( Pl + ir/2), 

P3 = ±(p 2 + vr/2), 
p 4 = ±(p3 + vr/2), 
Pl = ±(p 4 + vr/2), (31) 

where we ignore the arbitrariness of adding ±2n. As seen from these relations, the solution 
appears as an one-dimensional curve. For example, you find the following curve satisfies the 
above equations: 

(pi,P2,P3,Pi) = (p, ~2~ p ' P> ~2 ~^' ^ 

where p can take an arbitrary value. We can also understand this from the fact that (1271) 
reduces to three independent equations for a = 0. This emergence of a line solution indicates 
that the excitation around zeros for a small a becomes incompatible with a Lorentz-covariant 
or a Dirac form. Due to these nonlinear corrections, we conclude that we should use this 
action only with a substantially away from a ~ 0. This property at a = is also the case 
for any even number of dimensions. 

We can visualize the allowed range of the parameter a from the 2-dimensional analogy. 
In the 2d case the upper limit of a/2 is illustrated in Fig. HI For a = a/2 the two zeros join 
to a double zero, where the excitation is not Lorentz-covariant. For a lower limit, the 2d 
analogy shows "Minimal-doubling" persists as long as a is positive. We depict 2d cases for 
a = 0.1 and a = in Fig. [5]and|6j For a = the solution appears as an one-dimensional 
line, thus there is no simple pole in the propagator. 

We conclude that the parameter range for minimal-doubling in both 2 and 4 dimensions 
is given by 

< a < y/2. (33) 

In this range two zeros of the 4d Dirac operator, which we denote as p^ and p^ , are given 
by 

py = arcsin I J = arcsinA, (34) 



a-V2- 


-a 2 


2 


a + V^ 1 


- a 2 



~(2) ■ fa + V2 — c^\ . /- tk /0 ^ N 

py = arcsin I I = arcsin vl — A z , (35) 
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FIG. 4: A 2d analogy of two zeros at the upper limit a = y/2 in (|25p . Red-dotted and blue-solid lines 
stand for zeros of the coefficients of 71 and 72, namely (sin p\ + cos p 2 — \/2) and (sinp 2 +cospi — \/2) 
respectively. Two zeros collide with each other and reduce to a double zero shown as a black point, 
whose excitation is not Lorentz-covariant. 



where for convenience we define a parameter A as A = (a — a/2 — a 2 )/2. These two zeros 
reduce to Eq. ( fTBl with a = 1. Excitations from these zeros are given by 



D {1 \q) = iy/1 - A 2 (^ 1 q 1 + 7 2 g 2 + 73<?3 + 74^) - iA("/iq 2 + 72^3 + 73<?4 + 74?i) + 0(q 2 ), 

(36) 

£> (2) (g) = zA(7igi + 7 2 <?2 + 73<?3 + 74<?4) - Wl - A 2 (j 1 q 2 + 72^3 + 73<?4 + 74<?i) + 0(q 2 ), 

(37) 

where 0(q 2 ) corrections are not negligible for a small a as we have discussed. Momentum 



,W 



bases at the zero pjj, are given by b^ = (a/1 — A 2 , 0, 0, —A), b 



(i) 



J 1 ) 



>(!) 



(-A,VT^P,0,0), 

(0, -A, a/1 -A 2 , 0) and b^ 1} = (0, 0, -A, a/1 - A 2 ) while those at the other zero are 



.( 2 ) 



(2) 



(2) 



given by b^ = (A, 0, 0, -y/l - A 2 ), \%> = (-y/1 - A 2 , A, 0, 0), t# ; = (0, -Vl - A 2 , A,0) 
and b^ = (0, 0, —a/1 — A 2 , A). They are clearly non-orthogonal for cases of a 7^ 1 thus the 
associated lattices are non-orthogonal as with the Creutz fermion [l3| . 

We note we have confirmed the parameter range for minimal-doubling ( |33l) is also the 
case with other even dimensions (d = 2, 4, 6, 8, 10). From this fact we can anticipate this is 
a universal property for this class of actions for any even dimensions. 

Here we comment on the ranges of the parameter for minimal- doubling in odd dimensions. 



11 



a=0.1 
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FIG. 5: A 2d analogy of two zeros for the lower limit a = 0.1 in (|25p . Red-dotted and blue-solid lines 
stand for zeros of the coefficients of 71 and 72, namely (sinpi+cosp2 — 0.1) and (sinp2+cospi — 0.1) 
respectively. Two zeros shown as black points go further to each other with a decreasing. There are 
nonlinear corrections to the dispersion relation, which leads to Lorentz-non-covariant excitations. 
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FIG. 6: A 2d analogy of two zeros for the lower limit a = in (|25[) . Red-dotted and blue-solid 
lines stand for zeros of the coefficients of 71 and 72, namely (sinpi + COSP2) and (sinp2 + cospi) 
respectively. There are piled lines composed of red and blue lines, which stand for one-dimensional 
line solutions. 
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For example, in the 3d case, we can rewrite analogous equations for zeros of the operator 
with sinpx = x, sinp 2 = V an d sinp 3 = z as 

x 2 + y 2 — lax = 1 — a 2 , 
y 2 + z 2 — 2ay = 1 — a 2 , 



z + x — 2az = 1 — a . 



By eliminating y and z, we obtain the following equation for x, 



(38) 



(2x 2 - lax + a 2 - l)(8x 6 - 24ax 5 + (36a 2 - 12)x 4 - (36a 3 - 24a)x 3 

+ (2a 4 - 40a 2 + 6)x 2 + (10a 5 + 28a 3 - 6a)x + 25a 6 + 29a 4 + 11a 2 + 1) = 0. (39) 

As seen from this there are two real solutions similar with the 4d case for a < v2 as 



a ± a/2 - a 2 
x = y = z= . (40) 

However, in this case we have six more complex solutions which come from the other factor 
( 139]) . The question is how many real solutions this equation has. We have studied it 
numerically and found six real solutions exist below a critical value of a as a = 0.28... while 
there are six complex solutions above it. Thus the minimal-doubling parameter range in the 
3d case is given by 

0.28... < a < y/2. (41) 

In this case the lower limit is non-zero while the upper limit is the same as even dimensional 
cases. Indeed we can numerically confirm this kind of nonzero lower limit is also the case 
with any odd dimensional twisted-ordering action. We can explain the reason that there 
is such difference between the even and odd dimensional cases as following: As seen from 
(158"]) . there are d independent equations even for a = in the (d — 2n + l)-dim cases 
(n = 1, 2, 3, ...) while there are only d — 1 equations in the even dimensions for a = 0. This 
means that the Dirac operator can have 2 d zeros for a = in the odd dimensional cases 
while zeros appear as a curve in even dimensions. Actually, as seen from ( 139]) . there are 2 d 
independent real solutions for a = in the odd dimensions. Since only two zeros exist for 
a = 1 also for odd dimensions, there must be a critical point between a = and a = 1 in 
odd dimensions. 
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IV. SYMMETRIES 

In the previous section we have studied two classes of minimally doubled actions in ( 120]) 
and ( l25j) . both of which are obtained from the twisted-ordering action in (1121) . In this section 
we will discuss discrete symmetries of these actions and redundant operators generated by 
loop corrections. 

Firstly we note both of the actions possess some common properties with other minimally 
doubled fermions: "gamma-five Hermiticity" , "discrete translation invariance", "flavor- 
singlet £/(l)y" and "flavor-nonsinglet U(1)a" The latter is the exact chiral symmetry pre- 
venting additive mass renormalization for the neutral pion. In addition to them, they have 
exact locality and gauge invariance when gauged by link variables. We expect every sensible 
minimally doubled fermion to possess these basic properties. 

On the other hand, Discrete symmetries associated with permutation of the axes and 
C, P or T invariance depends on a class of the actions. In other words, we can identify 
minimally doubled actions by these discrete symmetries. Here we show discrete symmetries 
which the dropped twisted-ordering action possesses: 

1. CP, 

2. T, 

3. Z2 associated with two zeros. 

Here we take the p^ direction as time. We also write the symmetries which the full 
twisted-ordering action possesses: 

1. CPT, 

2. Z4 associated with axes permutation, 

3. Z 2 associated with two zeros. 

There is also the possibility to utilize the above Zi symmetry in order to recover flavored- 



C, P or T symmetry. Actually Borici-Creutz action has a flavored-C symmetry [17|. How- 
ever we do not further discuss here details of this point, to which a future work will be 
devoted. 

The dropped twisted-ordering action has CP and T invariance. We consider what redun- 
dant operators can be generated radiatively in this action following the prescription shown 
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in [171 ] . The relevant and marginal operators generated through loop corrections are given 

by 



Of = ^71^, 
of = ^72^, 
C>f = ^73^, (42) 



and 



Of = fty^, 
Of = ^iD^, 

Of = fy^D^, 

Of = faD^, 



of 


" \iv" \JLV1 


of 


= F lu F lu , 


of 


= F 2u F 2l/ , 


of 


= F2, v F- iu . 



(43) 

These operators are allowed by the symmetries for possessed by the action. Note that the 
marginal operators include a renormalization of the speed of light, both for the fermions and 
the gluons. 

On the other hand, the case of twisted-orderin g ac tion is similar to that of Borici-Creutz 



action which has higher symmetry S4 than Z4 17J. It is unlikely that there are more 



redundant operators generated in the twisted-ordering action than those of Borici-Creutz 
one shown in [l7J . However in the future we need to confirm that this reduction of symmetry 
does not affect the number of redundant operators. Then we will know whether or not this 
class is more useful than the Borici-Creutz class. 



V. CLASSIFICATION 

In this section we discuss a classification of the known minimally doubled fermions. So 
far we have seen four variations on minimally doubled fermions: Karsten-Wilczek, Borici- 
Creutz, twisted-ordering and dropped twisted-ordering fermions. Here we classify these 
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actions into two types: One of them, which includes Karsten-Wilczek and dropped twisted- 
ordering actions, is given by 



d (p) = X/ 7 ^ sin ^ + y^JliRjjjcosPj - 1] 



(44) 



i-.i 



where we take % = 1,2,3 as i while j = 2,3,4 as j. The point is that indices % and j are 
staggered. The different actions in this class depend on the choice for the matrix R. For 
example, consider the following ffs 



R 



R 



(l 


1 


l\ 











1° 





V 


(l 





°] 





1 





,0 





1, 



(45) 



(46) 



R 



(l 


i o\ 





1 


v° 


0, 



(47) 



In the case of (145]) . the general form (T44]) reduces to Karsten-Wilczek action as following, 

D{j>) = (sin pi + cosp 2 + cosp 3 + cosp 4 — 3) r/i 

+ sinp 2 «72 

+ sinp 3 z7 3 

+ sinp 4 z74. (48) 

Here you can also add an overall factor as a parameter to R. Actually the original Karsten- 
Wilczek action includes such a parameter A in front of 0(a) terms with a minimal-doubling 
range of the parameter A > 1/2. 

In the case of ( 146]) . this reduces to the dropped twisted-ordering action 

D(p) = (sinp! + cosp 2 _ 1) iji 
+ (smp 2 + cosps - 1)272 
+ (sinp 3 + cosp 4 - 1)273 
+ sinp 4 z7 4 , 



(49) 
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where you can also add an overall parameter to R. 

Finally, the action for the case (J4"T|) is a new possibility, given by 

D(p) = (sinpi + cosp 2 + cosp 3 — 2) i'ji 
+ (sinp 2 + cos£> 4 - 1) «7 2 
+ sinp 3 273 
+ sinp 4 z7 4 . (50) 

We see that there are many options associated with possible i?'s. One common property 
with fermions obtained in this way is that two zeros are separated along a single lattice axis 
and given by (0, 0, 0, 0) and (it, 0, 0, 0). One significance about the general form (I44J1 is that 
a coefficient of at least one Gamma matrix has no (cosp^ — 1) term. It is also notable that, 
in a coefficient of each gamma matrix, a momentum component p^ associated with a sine 
term differs from the component in any cosine term. These two points seem to be essential 
to minimal-doubling. Here we note minimal-doubling still persists if we add a (cospi — 1) 
term to any actions obtained from ( )44|) . We can change the location of the two zeros by 
adding such a term with a parameter. 

Here let us comment on a sufficient condition of R for minimal-doubling. We found 
minimal-doubling is realized when you take the following class of -R's: i?'s composed by 
three column vectors as R = (vi,v 2 ,v 3 ) with each vector Vi (i = 1,2,3) taken as (1,0, 0) T , 
(0, 1, 0) T or (0, 0, 1) T . All the examples we have shown here satisfy this condition. However 
we have not yet revealed a necessary condition for this. Our future work will be devoted to 
revealing such a condition. 

The second class of actions includes the twisted-ordering fermion with the a parameter 
and the Borici-Creutz actions. A generalized Dirac operator for this type is given by 

d (p) = 2^[7 M sin(p M + /3 / J - 7^ sin(p M - /3 M )] - iT, (51) 

A* 

where 7' = A^y is another set of gamma matrices where we define A as an orthogonal ma- 
trix with some conditions: At least one eigenvalue of A should be 1 and all four components 
of the associated eigenvector should have non-zero values. Here (3^ and T has a relation with 
this A as T = ^7^8^2/3^ = X]u7/i s i n 2/3^, which means sh^/^ is an eigenvector of A as 
A^y sin 2(3 U = sin 2/3^. Thus once A is fixed, /3 M and T are determined up to a overall factor 
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.4 



(52) 



of sin 2/3^. By imposing these conditions on A, /3 M and T, the action (IBTj) can be a minimally 
doubled action. In such a case /3 M indicates locations of two zeros as p M = ±/9 M . Adjusting 
/3, we can control the locations of the zeros. 

Note that the first and second terms in ( |5TT) are nothing but naive fermion actions. We 
can eliminate some species by combining two naive actions with different zeros in one action. 

Now let us show this action includes minimally doubled actions. Firstly this general form 
in (15TJ) reduces to Borici action by choosing A, which is given by 

7^ = A^u, 

/-l 1 1 1 \ 

1-111 

11-11 
\l 1 1-1/ 

where we fix (3^ = 7r/4 thus T is given by T = J7 sin(7r/2)7 M = ^«7m- 

This general form (15TJ) also yields the twisted-ordering action. For this we take the 
following matrix A 

7^ = A^ u , 

/o l\ 
10 
10 

\0 1 0/ 

fix Pn = it/A, and T is given by T = ^2 7m- Then the general form reduces to the twisted- 
ordering action with a = 1. We can also obtain the q^I case by choosing other values of 
/3 M . For < a < \^2 we take < (3^ < 7r/2 along with which an overall factor for T changes. 
We note this general form ( 15TJ) also includes non-minimally doubled actions like the 
partially twisted-ordering actions in ([8]) and ffTUl) . Here let us consider 



,4 



(53) 



7 



.4 



■^■flVIUI 

/o 1 o o\ 

10 

1 

\0 1 0/ 



(54) 
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and fix sin2/3 M where T is given by T = XL 7m- This clearly gives the partially-twisted- 
ordering action in ([TO]) . 

Finally we have classified the known minimally doubled fermions into the two types. 
Now we can derive a lot of varieties from these general forms deductively as we have shown 
in (150]) . However the second type of general actions (1511) is not restricted only to mini- 
mally doubled fermions. We may be able to restrict it by imposing other conditions. We 
note the two fermion actions which we derived from the twisted-ordering method, namely 
twisted-ordering and dropped twisted-ordering, belong to different types of minimally dou- 
bled fermions although they look similar to each other. 

VI. SUMMARY 

In this paper we study a variety of new classes of minimally doubled fermions and a 
classification of all the known cases. In Sec. Hi] we propose a systematic method to reduce 
the number of species, the "twisted-ordering method". In this method we can choose 2, 4, 
8 and 16 as the number of species by controlling discrete symmetry breaking. By using this 
method we obtain new classes of minimally doubled fermions, which we call twisted-ordering 
and dropped twisted-ordering actions. In Sec. HHI we study these two classes in terms of the 
parameter range for minimal-doubling. In Sec. IIVI we study discrete symmetries of these 
fermion actions and show twisted-ordering action has Z4 symmetry while dropped twisted- 
ordering action possesses CP invariance. We also show they require fine-tuning of several 
parameters for a correct continuum limit because of lack of sufficient discrete symmetries in 
the known classes. In Sec. |V] we classify all the known minimally doubled actions into two 
types. We can derive several unknown minimally doubled actions by this classification. 

With several varieties of minimally doubled fermion actions available, a goal is to apply 
these actions to numerical simulations and study their relative advantages. In particular 
it is important to investigate how many fine-tuning parameters are required for a good 
continuum limit and how difficult this tuning is. This kind of study was started in [17]. 
Recently renormalization of the operators has been also studied in 22j-|24j. Both of these 
studies focused on Karsten-Wilczek and Borici-Creutz fermions. One can now explore other 
varieties such as ones in this paper. With a better understanding of their properties, one 
can go on to larger numerical calculations. 



22 



Acknowledgments 

TM is supported by Grand-in-Aid for the Japan Society for Promotion of Science ( JSPS) 
Research Fellows(No. 21-1226). TM thanks Taro Kimura for stimulating discussions. MC is 
grateful to the Alexander von Humboldt Foundation for support for visits to the University 
of Mainz. This manuscript has been authored under contract number DE-AC02-98CH10886 
with the U.S. Department of Energy. Accordingly, the U.S. Government retains a non- 
exclusive, royalty-free license to publish or reproduce the published form of this contribution, 
or allow others to do so, for U.S. Government purposes. 



[1] K. G. Wilson, Phys.Rev.D 10, 2445 (1974). 

[2] H. B. Nielsen and M. Ninomiya, Nucl. Phys. B 185, 20 (1981); Nucl. Phys. B 193 173 (1981); 
Phys. Lett. B 105 219 (1981). 

[3] D. B. Kaplan, Phys. Lett. B 288, 342 (1992). 

[4] V. Furuman and Y. Shamir, Nucl. Phys. B 439, 54 (1995). 

[5] P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25, 2649 (1982). 

[6] N. Neuberger, Phys. Lett. B 427, 353 (1998). 

[7] L. H. Karsten and J. Smit, Nucl. Phys. B 183, 103 (1981). 

[8] J. B. Kogut and L. Susskind, Phys. Rev. D 11, 395 (1975). 

[9] L. Susskind, Phys. Rev. D 16, 3031 (1977). 

[10] H. S. Sharatchandra, H. J. Tlnm and P. Weisz, Nucl. Phys. B 192, 205 (1981). 

[11] L. H. Karsten, Phys. Lett. B 104, 315 (1981). 

[12] F. Wilczek, Phys. Rev. Lett. 59, 2397 (1987). 

[13] M. Creutz, JHEP 0804, 017 (2008) [arXivi0112j20TJ [hep-lat]]. 

[14] A. Borici, Phys. Rev. D 78, 074504 (2008) [arXiv:0712.440Tl [hep-lat]]: PoS LATTICE2008, 

(2008) [arXiv:0812.0092] . 

[15] M. Creutz, PoS LATTICE2008, (2008) |arXiv:0808.00T4] . 

[16] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim, Rev. Mod. Phys. 

81, 109 f2009h larXTv:0709.1163l [cond-mat.other]. 

[17] P. F. Bedaque, M. I. Buchoff, B. C. Tiburzi and A. Walker-Loud, Phys. Lett. B 662, 449 (2008) 



23 

|arXiv:0801.336H [hep-lat]]; M. I. Buchoff, PoS LATTICE2008 (2008) jarXiv : 0809. 3943 j . 
[18] P. F. Bedaque, M. I. Buchoff, B. C. Tiburzi and A. Walker-Loud, Phys. Rev. D 78, 017502 

(2008) [arXiv:0804.fT45l [hep-lat]]. 
[19] T. Kimura and T. Misumi, to appear in Prog. Theor. Phys. [arXiv:0907.1371 [hep-lat]]. 
[20] T. Kimura and T. Misumi, Prog. Theor. Phys. 123, 63 (2010) [arXiv:0907.3774l [hep-lat]]. 
[21] K. Cichy, J. Gonzalez Lopez, K. Jansen, A. Kujawa and A. Shindler, Nucl. Phys. B 800, 94 

(2008) [arXiv:0802.3"637l [hep-lat]]. 
[22] S. Capitani, J. Weber, H. Wittig, Phys. Lett. B 681, 105 (2009) (arXiv:0907.2825j . 



[23] S. Capitani, J. Weber, H. Wittig, (2009) [arXiv:0910.2597j . 

[24] S. Capitani, M. Creutz, J. Weber, H. Wittig, (2010) [arXiv: 1006.2009] . 



